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Abstract 

We consider two types of the time-dependent Ginzburg-Landau equation in 2D bounded 
domains: the heat-flow equation and the Schrodinger equation. We study the asymptotic 
behaviour of the vortex solutions of these equations when the vortex core size is much smaller 
than the inter-vortex distance. 

Using the method of the asymptotic expansion near the vortices, we obtain the systems 
of ordinary differential equations (ODEs) governing the evolution of the vortices. The ex¬ 
pressions for these equations in the circle and in the annular domain are presented. We study 
the motion of the vortices in these two domains. It is shown that there exist the stationary 
points for both types of the equations and these points are determined for some particular 
cases. 

The heat-flow equation describes the particles that move like electric charges. If the 
vortices have different signs (i.e. Poincare indices) then they attract each other and collide. 
If they have the same signs then they repulse. The particles always move away from the 
stationary points. The motion is very slow and non-periodic. 

The Schrodinger equation describes the motion of particles that behave like hydrody¬ 
namics vortices. The vortices with the same signs move in the same direction. If the signs 
are different then the vortices move in the opposite directions. In particular, if the initial 
positions are near the stationary points, then the particles move along the elliptic trajec¬ 
tories. The motion is not stable with respect to the initial data. It is always periodic or 
quasiperiodic. Examples of such trajectories are presented. 
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1 Introduction 


We start with the free energy in Ginzburg-Landau theory: 

F {u) = + a\u\ 2 + ^\u\ A . ( 1 . 1 ) 

Here in 0) u{x.,t) (x = (x,y)) is a complex-valued wave function, H is the domain occupied 
by (super)fluid, m is an atomic mass, 2-kH. is a Planck’s constant, a and (5 are the temperature- 
dependent parameters such that 


a(T) = chq 



(3{T) >0, a 0 > 0 


with T c being a temperature of phase transition, i.e. a < 0 for T < T c . 

In the absence of gradients F(u) represents the condensation energy with the density 
f(u) = ot\u\ 2 + t)/3|u| 4 . 

We assume that there are N vortices in the domain, located at the points £i(i),... ,£jv(i)> = 

(Q, fjj) with degrees ni,... ,njv, respectively. The vortices are stable if and only if rij = ±1. If e 
is a size of vortex core, than inter-vortex distance is supposed to be of the order 0(1). 
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Let us pass to dimensionless variables u = uqu', x = ex. 1 , where ttg = — — corresponds to the 

r' 


minimum of the condensation energy and e = 
core). Then the free energy takes the form 


K- 2 


2m|o;| 


is a coherence length (size of vortex 


F(u')= [ |W| 2 + 

Jn 2 


( 1 . 2 ) 


Note that in these variables the vortex core size is 0(1) and inter-vortex distance is 0(l/e). 

Using standard techniques of the calculus of variation, the minimization of F with respect to 
variation in v! yields the stationary G-L equation: 


A u' + u\ 1 — \u'\ 2 ) = 0. 


(1.3) 


In this paper we study two types of non-stationary Ginzburg-Landau equations: the nonlinear 
Schrodinger equation 

A v! + u' {1 — | u'\ 2 ) = —iu[] 


and the nonlinear heat flow equation: 

A u' + u' {1 — \u'\ 2 ) = u' t . 


(1.4) 

(1.5) 


Under the usual diffusive scaling x" = x e, t" = e 2 t, Eq. (1.4) takes the form 


A u e + —u £ ( 1 — | |") = —iul 


and for (1.5) we obtain 


Au £ + A;u £ ( 1 - |rt e | 2 ) = u\ 


( 1 . 6 ) 

(1.7) 


We have again the vortex core size 0{e) and inter-vortex distance 0(1). 

In this paper we consider these equations in the bounded domain with the Neumann boundary 
condition on <9U. 

We look for two approximate solutions of ( |1.6|) and O) (see 0): outer solution 

v £ (x, t ) = rto(x, t) + etti(x, t) + e 2 rt 2 (x, t) + ... (1.8) 

in the region outside the vortex cores, but for small r = |x — £j(t )|, and inner solution for the 
core of the j th vortex: 


E7‘(x, t) = Uo + sUl ( + ^ ( iAifi) , + ... 


(1.9) 


Substituting the expansions (1.8) and ( |L9| ) in the equations (0) and (Q) , we obtain the 
series of equations for the functions ito(x,f), rq(x,t), 7x2(x, t),..., and Uo(x,t), U\(x, t),.... 
We construct the approximate form for the solution 0 ) near the jth vortex. For the inner 

solution (0) we pass to the stretched variable X =- anc [ the solution for large 

R = |X|. Asymptotic matching of these solutions yields the equations governing the motion of 
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the particles. Taking into account that r = eR, we obtain that the vortices described by the 
Schrodinger equation are gouverned by the system 

ij(t) = 2KW, (1.10) 

where the vector depends on the form of the domain O and can be expressed, for example, 
through the derivative of the complex potential W(z) of the fluid with vortices out of jth vortex: 


K (j) = (ReW'(zj), —ImW'(zj)). 

In particular, for annular domain we have (see formula (|2.10|) below): 


( 1 . 11 ) 


N 


W'izj) = — J2 n k 


k=1 


Zk 


C(*ln^)- C f*ln^ 




-^C(2iln-^) + 


2 ir\ 


N 


"■‘it: 




We obtain by the same matching procedure for the heat-flow equation that for the standard time 
scaling t" = e 2 t the vortices are stationary at the leading order. To obtain a non-trivial dynamics 
we have to rescale the time by the factor (logl/e). In this long-time scaling the dynamics of 
the vortices is described by the system 

if{t) = 2n^\ ( 1 . 12 ) 


where = (— 7 ) j,Cj) and K^ is determined by ( 1.11 ). 

We construct the numerical solutions of the systems (1.1C) and (1-12) in the circle and in the 
ring for given initial positions of vortices and given indexes. We show that the Schrodinger 
equation (1.10) discribes the motion of the hydro dynamical vortices and the heat-flow equation 
( 1.12 ) governes the motion of some charged-similar particles. The trajectories of the particles 
for different initial data are shown in the Figs. 1-9. 


2 Schrodinger equation 


Let us consider the Schrodinger equation (L6). Substituting (L8) into (0) , we obtain the 
following the series of equations: 


e 2 : u 0 (l - K| 2 ) = 0; 


( 2 . 1 ) 


e 1 : u\Uq + uqU\ = 0; 

e° : — iuot = Auo — |ui| 2 uo. 

Here u means the complex conjugation of u. The first of these equations yields 

uo(x,t) = e i ®o( x ,t). 

Now from (2.3) and (2.1) we get 

•Lot = iAd > 0 - (V<L 0 ) 2 - \ui\ 2 . 


( 2 . 2 ) 

(2.3) 

(2.4) 
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Separating real and imaginary part, we obtain 


Ad> 0 = 0. 


Asymptotic boundary conditions at the cores are 


d>o(x, t) —■> rij tan 1 j = nj6 

as x —> £j, where (r, 9) are the polar coordinates related to the jth vortex: r 

-■-(Kf 

Hence, a boundary-value problem for &o(x,t) takes the form: 




A<f> 0 (x,t) = 0, 

$ 0 CM) —> rijO, 

d&o 


in n \ 


x —> j = 1,2,... ,N 


dn 


= 0 


x 


G dVt. 


(2.5) 


The last boundary condition follows from the condition for u £ on <9fh It is well-known (see e.g. 
Q) that the phase of condensate wave function <f>(x, t) is related the hydrodynamic velocity 
potentilal ^>(x, t) 

0( x ,t) = — <f>(x,f). 


m 


h. 


Then we obtain v = = —Vd>, and the boundary condition is the standard hydrodynam- 

m 


ics condition of impermeability, v • n = 0. Thus we can consider problem (2.5) as a usual 
hydrodynamic problem describing the potential of the ideal liquid with the vortices. 

The solution of this problem can be found for any form of the domain H by the method of 
“reflexed vortices” or by some other method. In particular, for the annular domainith < r < R 2 
the velocity potential <ho( x , t) of an ideal fluid with N vortices of vortex strengths rik can be 
written in the complex form (see. ||) 

( N 


3>o(z) = ReW ( z ) = Re < — i nk 

l k=1 


In a [ i In — 

Zk 


— In a 



2 r?, r k 
— In —— In z 

LO\ R2 


+ const. 


( 2 . 6 ) 


Here z is a complex variable: z = x + iy = r exp(i6), z k = Cfc + = r k exp (i0 k ) are the points 

of vortex locations, r = |x|, the polar angle 0 = tan -1 (^) is measured from the center of the 
ring, W(z ) is a complex potential, Z k = R 2 exp (i0 k ), ((z) and a(z) are the Weirstrass zeta- and 
sigma-function with the half-periods ui\ = ir, u >2 = iln(R 2 /R\); r) = ^(cji).(see e.g. [§]). 

To obtain the approximate expression for the wave function uq near the j th vortex, let us expand 
<h 0 in a Taylor series around jth vortex. We have to take into account that $0 has singularities 
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rij tan 1 (y — rjj)/(x — Cj) at the points of vortex locations £j. Extracting this singularitity, i.e. 
representing <ho( x i t) = rijO + 'Loj ( x ; t) and expanding <3?oj( x > t) in the series 

$0j( x >*) = ®oj(£j,t) + V$ 0 j(£j,t) • (x-^) + ..., (2.7) 


( 2 . 8 ) 


we obtain the relation 

u 0 (x, t) = e^° « e in i e+ie °[l + iK® • (x - £,) + 0(r 2 )], 

where 

Oo = ®0j(£j,t) 

and 

K® = V$ 0 = (ReW'(zj), —ImW'(zj)). (2.9) 

Here W'(zj) is a derivative of the complex potential in the point of the vortex location. 

In particular, for annular domain we have (see [fl) 


1 N 

W\ Zj ) = - 


k =i 


Zk 


C iIn — ) — £ (i In Z ^ Zk 


zl 


N 


—C ( 2iln ih ) +—^S^ln^r- (2.10) 


One can show that the expression ( |2.10 ) gives the velocity potential of the ideal liquid with the 
vortices it the circle when R\ —* 0 (see Section 5) 


N 


»"(%■) = E 7 


nk 


N 


E n k 
it *_• — 


i{zj ~ Zk) ^ i(zj - 4) 


( 2 . 11 ) 


i? 2 

where z' h = —- are the coordinates of the vortices symmetric to Zk with respect to the circum- 

Zk 

ference r = i? 2 - 


Thus, the exteriour solution uq( x T) is defined by ( |2.8D , (2.9) and (2 .10| ) (or (2 . ll|) ) 


Let us consider now the core structure of the vortices. We denote by X =- — - - the stretched 

variable and by (R, d) are the polar coordinates of X. Then the inner solution^ Q takes a form 

*7 £ (x,f) = t/ 0 (X) + el7i(X) + e 2 t/ 2 (X) + ... 


Substituting this expansion in the equation (1.7) we obtain the series of equations for t/o(X), 
t/i(X), t/ 2 (X)...: 

AU 0 + U 0 (1-\U 0 \ 2 ) = 0, (2.12) 


A U x + (1 - 2\Uq\ 2 )U\ - UqU\ = i£j(t ) • W 0 , (2.13) 

As in the consideration above (see. ( |2.8D ) we look for the solution of ( |2.12| ) in the form: 

U 0 (R, 9, t) = f 0 (R)e in i d+ie °. 


Then fo(R) satisfies the equation 

1 Tl 2 

/o+^/o-^/o + /° (1 -/o 2) = 0 


6 






















with boundary conditions 

f/o(0)=0, 

\/o(oo) = 1. 

As in ||] it can be shown that 

1 -^ + °(si) (2 ' 14) 

for i? » 1. 

As for U i (X) = Ui(R, 6), we look for it in the form 
U X {R,6) = h{R,0)e in i e+ie ° = [(A r (R) + iAi(R)) cos 9 + (B r (R) + iBi(R)) smd]e in > 9+i9 °. (2.15) 


Substituting Uq and U\ in ( 2.13 ), gathering the terms with sin 0 and cos 6 and separating real 
and imaginary parts, we obtain the system of equations for A r , A ^, B r and B, : 

■ fo ah , 1 At , (-, o/2 2 \ 2 Bi 

- Ujr]j -R ~ Ar + ~R Ar + ( 1 ~~ 3 ^° _ ) Ar ~ 


- B r + ~ 3 /o 2 - -^2^ B r + n j-£2~i 


foCj — A” + —Af 


f/j. _ i ^ d/ i I 1 /2 


2 > 

| A, 

2 B r 

R 2 j 

+ & 

2 N 

)Bi 

2 A r 


n ’R 2 


.UO * R 1 

We look for the solutions A r , A t , B r , and B t in the form 

A r = -nji)jW(R), At = C 'jZ(R), 


B r = rijCjW(R), Bi = rjjZ(R). 


(2.16) 


Substituting these expressions into equations for A r -B{ we derive two equations for W(R) and 
Z(R): 


k = W " + k w > + w ( i - 3 / 0 2 - A) + A z , 


R 


R 


R 2 R 2 ' 


(2.17) 


f ° = Z " +l R Z , + W { l - f 2 °-l) + l? W - 


Taking into account ( 2.14 ), we conclude that 

f « = w +0 {i) 

for large R. 
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Expanding W(R) and Z(R) in the series with respect to the degrees of R, we can find the 
solution Ui(R,6,t ) for large R: 


Th(R,6,t) = e in > 


_ irijO+iQo 


rijC\ 

R 


X 


iR 


3 R ) + 2 


^ 3 R + R 




i 


where = {Cj,Vj)i £/ = i~VjXj)j Ci> -Ai are constants. 

Then the solution U (X) = C/q(X) + ef/j (X) can be represented as 


C/(X) = e' 


p irijO+i6o 


IF 

1 + -(^. X ) + o 


1 


1 + e e 


X 


/? 


X 

~R 
(2.18) 


(2.19) 


for large i?. 

The function [/(X) is an approximate solution in the region near the jth vortex, where R = 

l( x -£j)l/ e is lar S e - 

To guarantie consistency of the solutions (T8) and ( |2.1£ ) it is necessary to suppose that the 
leading terms of the expansion be equal for x — £ -(i) = cX, i.e. 


1 . 


• eX = -ij • eX. 


This gives the equations 


& = 2KW j = 1,2,..., N. 


( 2 . 20 ) 


This is the desired system describing the motion of the vortices in the domain fh 


3 Heat flow equation 


Now we consider the heat-flow equation O)- To construct the solution near the jth vortex in 
outer region we substitute (L8) in ([O]) and obtain the series of equations for uq , u\,. ..: 


u 0 (l - K| 2 ) = 0; 


uiuo + = 0; 

UQ t = Art 0 - |ui| 2 u 0 - 


For the function uq we have 

u 0 (x,t)=e i ^ t \ (3.1) 

as it was in the case of the the Scrodinger equation, but now <f>o satisfies the equation: 

$ 0i = A<f> 0 . (3.2) 

Passing to the moving coordinates x = x — fj(t) we obtain the following equation for <l>o that 
describes the local behaviour of the solution near the jth vortex: 

= A<h 0 (x,t), 

< (3.3) 

<ho(x, t) —► nj tan -1 ( y/x ) = Uj6 , x —> 0; 
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We denote by (r, 0) the polar coordinates corresponding to x. The function d>o near the vortex 
is described by the power series of r and logr for small r. It can be shown that 

<3?o(x,f) = njO + ^nj(\ogr)if -x + Sj • x + 0(r 2 logr), (3.4) 

where S j is an unknown vector. For the function no we have 

u 0 (x, t) = e njd [l + ^n,(log r)& -x + iSj • x + 0(r 2 logr)]. (3.5) 

z J 


This gives a representation of the outer solution near the jth vortex. 

The inner solution for the heat flow equation in the stretched variables leads to the system of 
equations for Uq, U \,..., similar to (2.12)-(2.13): 


(2.12)-( 2.13 


AU 0 + U 0 (1-\U 0 \ 2 ) = 0, 


(3.6) 


• W 0 = A lh + (1 - 2\Uq\ 2 )U\ - U^IR. 


(3.7) 


The hrst equation is the same as ( 2. 12|) a nd hence we have again (2.14). We look for the function 
t/i(X) in the form (|2.15|) . Then Eq.(|3.7|) in terms of /o, /i becomes 


-ij • {f 0 VR + injfoVO) = A h + 2in j (Vf 1 ■ W) - -^/i + /i(l - 2/ 0 2 ) - f 2 f v 


(3.8) 


The following system for A r (R), Ai(R), B r (R), Bi{R ) is valid: 


~Cjfo ~ A” + —A' r + fl - 3 /q - A r - nj-jA, 


1 


~Vjfo - B” + —B' r + ( 1 - 3 fo - -R 2 ) Br + n j-^Ti 


R 2 


2 Ai 


fo 


~ n j( jl R ~ + ( 1 - fo ~ tp> ) Ai + n 


10 h) 


R 2 
2 B r 

2 A r 


fo n 1 / ( o 2 \ 2^4 r 

njVj R ~ Bl+ Ti> + { 1 ~ fo ~lp) Bi ~ nj lP 


R 2 ' 

We seek the solution of this system in the form that is similar to 

A r = (jZ, Ai = fjjUjW ., B r = r/j Z. Bi = —QrijW, 
where Z and W satisfy 


~fo = Z" + j^Z 1 + ( 1 — 3/q — 
~R = W " + ^ W ' + [ 1 -fo- 


2 \ 2 W 

Z + 


R 2 

2 


R 2 ’ 

2Z 


, W + 

R 2 R 2 
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The solutions W and Z can be found in the power series of R and logit!. We are interesting in 
the behaviour of the solutions for R S> 1: 


W = —-RlogR + CqR + 0(log R), 


Z = -—logR + 0 


Here Co is a constant. So, we have for R 1 

1 


tfi(X) = e* 


^irij(log R){£j ■ X) + injCodj ■ X) + O 


logi? 

R 


(3.9) 


Putting together the solutions Uq and U\ we obtain as R —> oo 


in-jO 


U 0 (X) + eU 1 (X) = e m i 


fo(R) + (log R)(if ■ X) + ierijCodf ■ X) + O ( e 


log R 

~r~ 


(3.10) 

Matching ( |3 .1 0| ) with the outer solution (3.5) and taking into account that x = eX, r = eR, we 
can see that 

^(loge)^ + rijSj = Coif 


or 


Zj = —r 


2 v—o -/v? 

S ; 




-4 log e + C 0 


= 0( l/log(l/e)). 


(3.11) 


This equation shows that in this time scale the vortices are stationary at the leading order. The 
time scale of the vortices is slower than the time scale of the phase field. 


4 Heat flow equation: long time scale 


To obtain a non-trivial dynamics of the vortices, we have to rescale the time variable by a factor 
0 (log(l/e)), i.e. instead of diffusive scaling x" = x e, t" = e 2 t we consider the scaling x" = x e, 
t" = [e 2 /log(l /s)]t. Then we obtain the equation (see pj]) 


A u £ + -^u £ ( 1 - 


u £ \ 2 ) = 


log(l/e) 


Ut 


(4.1) 


We consider the approximate solutions of this equation in the form (1.8) for outer region and 
(l.£) for inner region. Thus we obtain as in the previous case the equations ( j2.4| ) for uq. For 
the function <f>o we have 

-$o t = A$ 0 - (4.2) 


1 


log(l/e) 


To investigate the behaviour of the outer solution near the j th vortex we expand Tn in powers 
of 5 = (log(l/e)) _1 : 

4*o = 4>oo(x,t) + (5$ 0 i(x,t) + (5 2 4> 0 2(x,f) + • • • 


Then we get the series of equations for 4> 0 j: 

A4>oo = 0, 


(4.3) 
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(4.4) 


$00 ,t = A$oi, • • • 

Eq.([hj|) shoud be completed by the boundary conditions at the vortices: 


<h 0 o(x, t) —> rij tan 


-l ( y~Vj 

x - Cj 


and the Neumann condition on the boundary i.e. 


<9$ 


oo 


dn 


= 0. 


Thus, near the jth 


So, we have problem (2.5) for $oo( x > t) and hence $oo( x ) t) is given by (|_ 

vortex 

$oo (M) = rijd + 0 o + K (j) • (x - 
where K® is defined by (2 .1C| ). 

For the function $qi we have 

$oi CM) = 7 ,™j (log r)(£j- ■ x) + P i j • x + 0(r 2 logr) 

for r = |x| < 1, Pij is a constant vector depending on the history of the vortices. Extracting 
the singularitity of < I > oi(x, t) in the jth vortex (see @3)) , we obtain 


uq( x , t) = e mjd ( 1 + iK^ ■ x + 


in -> (log r)(£j- ■ x) + (logr)(Pi j • x) ) + 


21og(l/e) 3 log(l/e) 

+0(5 2 + 5r 2 log r + r 2 ) 


(4.5) 


as r is sufficiently small. 


It follows from (4.1) that for the inner solution 

6 co 


log(l/e) 


VU 0 = AUi + (1 - 2\Uq\ 2 )Ui - UqU\. 


By the same way as it was done for (T7) we obtain 

U 0 (X)+eUi(X) = e in i° \f 0 (R ) + - . . (log R)(& • X) + ■ X) + O 


21og(l/e) 


log(l/e) 


as R —> oo. 


Matching (|T5|) and (|47| ) gives 

K (i) + 


-(log e)£j- + 


m 


log(l/er) 


= njCo - ± 
lj log(l /ey 3 ' 


(4.6) 


f glog R \ 
\R\og(l/£)J 
(4.7) 


(4.8) 


21og(l/e) 

Therefore we obtain the desired equation 

4 ± = 2n,K^ + 0(l/log(l/e)). 

In this time scale the dynamics of vortices is not trivial and it is described by the system of 
equations 

if = 2n i K«), j = 1,2,..., N . (4.9) 

This equation is similar to ( 2.20| ) for the Schrodinger equation, but it includes the turned co¬ 
ordinates and the factor n.j. Notice that the time scale in (4.9) is slower by the factor 
0( log(l/e)). 
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5 Motion of the vortices in the circle 


Let us show that from (2.6) we obtain the velocity of the liquid with the vortices in the circle: 


N 


W'(z) = £ - 

Z—/ 7 


nk 


N 

S77 




^ - z k) K z - 4) ’ 


(5.1) 


R 2 

as R\ —> 0. Here 4 = ~ are the coordinates of the vortices symmetric to Zk with respect to 

Zk 

the circumference r = Rq. 


It follows from (2.6) that 


N 


w l (z) = -T 


nk 


k =1 


Zk 


C iln — -C iln 


zz k 

zl 


-4( ak ^) + SE"*i» 


r \ 2 ir] 
i?2 J W\Z 


N 


k =1 


/L 


(5.2) 


Notice that as > 0 the denominator g tends to zero. Therefore, the ((-function 


. , N 77 £ 7T 71 ^ . v—> 

^) = - + — cot — + 4^ 


LO 1 2wi 


00 2 1 

qZt 


~q 


2t 


Sill 


7 TtZ 


is reduced in the sum of the linear part and cotangents. Hence 

N 


k =1 




2 2* 


cot - In — — cot - In 


* , ZZ k 


2 Zl 


Calculating all the cotangents, we obtain the expression 


W'(z) 




-g + Zk 
z(z - Zk) 


ZZ k + 

z(zzk - Zl) 


Denoting by Zk = Z^/zk the coordinates of the vortex symmetric with the vortex Zk with 
respect to the circumference r = R -2 (Zk = /? 2 exp( 0 fc)) and expanding the terms into the 
simplest fractions we derive the desired formula (|5.l|). 

To obtain the velocity of the jth vortex we have to substract the term corresponding to this 
vortex and then put z = Zj. We have 


N 


= E 7 


nk 


N 


E n k 
i( — ?!. 


Tf *(*i “ z k) “ i(zj ~ 4) 


(5.3) 




Let us investigate the motion of the vortices in the circle. The systems (2.2C) and (4.9) describing 
this motion can be represented as 


Cj ~ 2Kj\, 

Vj = ~ 2K j2, 7 = 1,2,...,IV 


(5.4) 
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for the Scrodinger equation and as 

iCj = ~^ n j^j2, 5 n 

[Vj = ~ 2 n j K ji, j = l,2, ...,N 

for the heat-flow equation. 

Here K j = (Kji,Kj2) = ^RW'(zj), 3fW'(zj)) and W'(zj) is defined by (|5 . 3|) . 

5.1 Schrodinger equation 

Notice that the particles described by the Schrodinger equation behave as hydrodynamic vortices. 
So, we can expect that they move along the closed trajectories and never meet the walls. 

If at the initial moment we have the symmetric locations of the vortices with alternating signs, 
we obtain the symmetric closed trajectories (see. Fig.l for N = 2 and N = 4). The closer are 
vortices to the wall at t = 0 the wider are trajectories. 


rl =(0.5,0.0), r2=(-0.5,0.0), n1=-n2=1 


r1=-r3=(0.8,0.4), r2=-r4=(-8.0,0.4), n1=-n2=n3=-n4=1, T=2.0 




Fig.la. Motion “vortex-antivortex”. Fig.lb. Motion ”vortices-antivortices”. 

N = 2. N = 4. 

If one of the vortices is close to the center and the other is near the wall we have more complicated 
motion: the vortex near the wall moves almost along the circumference and the inner vortex 
passes a “star” trajectory (Fig.2a). But if all the vortices have the same signs and are on the 
same circumference at t = 0, they move independently along this circumference (Fig.2b). 

The motion is much more complicated for odd number of the vortices with the different signes 
(Fig.3a, 3b): vortices move in the chaotic way even if the initial positions are on the same 
circumference (comp, with Fig.2b). 


5.2 Heat flow equation 

The figures below represent the motion of the particles governed by the heat flow equation. In 
fact, we have only two types of motions: attraction, if the vortices have different signs and in 
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rl =(0.1,0.1), r2=(-0.1,1.3), n1=-n2=1, T=8.0 


ri=(cos(2*p*i/N),sin(2*p*i/N)), ni=1 fo all i 




Fig.2a. Motion “vortex-antivortex”. Fig.2b. Motion of the same vortices on 

Non-symmetric initialy positions. the same circumference. N = 5. 


r1=(1.0,0.0), r2=(-1.0,0.0), r3=(0.0,1.0), n1=n2=-n3=1, T=3.0 


ri=(cos(2*pi*i/N),sin(2*pi*i/N)), n(2*k-1)=1, n(2*k)=-1, T=1 




Fig.3a. Motion “vortices-antivortices”. 
Symmetric initialy positions. 


Fig.3b. N = 5. 


the initial time are close enough to each other, and repulsion, in all other cases. 

So, if we have the even number of particles with alternative signs, each particle is attracted to 
the nearest neighbour and annihilated (Fig.4a, 4b). 

For the odd number of particles the picture is symmetric if the initial data is symmetric. In 
Fig.5a we have the circle of particles with alternative signs. In this case the nearest particles 
with different signs are attracted and annihilated. 

Finally, any number of particles with the same signs is attracted to the boundary and annihilated 
with the imaginary vortices (Fig. 5b). 

Thus the motion of the particles in the framework of the heat flow equation is always attractive 
or repulsive. There are no periodic and quasiperiodic trajectories. 
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ni=1 for even i and -1 for odd i 


ni=1 for even i and -1 for odd i 




Fig.4a. Motion of different-sign particles: 
symmetric initial positions 


Fig.4b. Attraction to the nearest 
neihbour. 


ni=1 for even i and -1 for odd i, N=5 


ri=0.6(cos(2*PI‘i/N),sin(2‘PI*i/N)), ni=1 for all i 




Fig. 5a. Motion of odd number of 
different-sign particles. 


Fig.5b. Attraction to the boundary. 


6 Motion of the particles in the ring 


6.1 Stationary points 


We investigate the motion of the vortices in the ring. The motion is governed by (5.4) and ( p! 
but W'(zj) is defined by (|2.1C|) . 

Consider the expression for W'(zj). 


A. Let us suppose that there are only two vortices: N = 2. Representing Zj = r ] e i0:i we obtain 


W'(zj) = 


Tk 


nkC ( i In — - {Oj - 9 k ) ) - n k ( ( i In - (6b - Q k ) ) - 


Rl 
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( 6 . 1 ) 


~ n 3 C 



2*7/ / ri r 2 

H-ni In ——|- n 2 m — 


U)l 


R2 


R2 



Here j, k = 1,2, k ^ j (if j = 1, then k = 2 and vice versa, if j = 2, then k = 1). 

Notice that uq = it and cu 2 = *ln(i? 2 /i?i) are the half-periods of the ((-functions. 

If we put n = r 2 = ro (vortices are on the same circumference) and take into account that 
((—z) = — C(z))> we obtain 


W\ Zj ) = - 


Zi 


nkC {Ok - 9j) ~ n k C ^2 i In - ( 6 j - 9 k )^j - nj ( [ 2,1 In + “( n i + n z) ln 


Using the addition formula for the ((-function (Q, p.635) we obtain 


r 0 


R 2 


C 2iln-f-(0j-0 fc ) =C 2iln-f + C (0 fc - %) +- 


^0 


«2 


1 


where T^z) is Weierstrass’ P-function (||). Thus 


W'{ Zj ) = - 


~{n k + rij) 




2 


V (2iln-g) -V (9 k ~9j) 
V(2iln^-) -V(9 k -9 j ) 


V (9k-9j) 

V[2i\n^\ -V(9 k -9 j ) 


Notice that the function V'(±LUi) = 0 vanicshes at the points of the half-periods ±(Uj, i = 1,2. 
Thus, if we put 


r 0 R 2 

2i ln — = — u; 2 = — i ln — and 0 k — 9j = uj± = ir the last term will be equal zero. Hence, for 

R 2 _ Ri 

r\ = r 2 = y/R 1 R 2 , 9j — 9 k = tt we have 


W'(zj) 


ni + ri2 , rjLU2 

- T) - 

Zj l LO\ 


n\+ n 2 iri 
Zj 2uq 



n\ = —n 2 , 
n\ = n 2 = 1. 


( 6 . 2 ) 


Here we have used the Legendre’s relation ryo; 2 — r{vj\ = 7tz/2; rf = C(cu 2 ) (||). 

Expression ([h^) means that there exist the stationary points r\ = r 2 = sjR1R2 , #1 — 0 2 = n 
for both equations if ni = — n 2 . Without loss of generality we can suppose that n\ = —n 2 = 1 
and that the stationary points are 

(x° 1 ,y° 1 ) = (y/Rjh, 0), (x 0 2 ,y%) = (-y/RJh, 0). (6.3) 


If we take the initial positions of the vortices at these points, the vortices will not move at all. 


B. Let us obtain the similar points for any even number of vortices 2 N. It was shown (||]) 
that if we have N vortices with the equal intensities nj = n located uniformly on the same 
circumference r = ro, i.e. Zj = r exp(i.9j), 9j = 2nj/N, j = 1, 2,... , N, the expression for W'(z) 
take the form 


W'(z) = 


n-j 


C * ln — — C N In 


ZZj 

Rl 


2 in r 

-ln — 

UJ\ R2 


(6.4) 
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where Zj = roex.p(i9j ) is the position of any vortex and the half-periods of the (-function are 
cui = n/N and ui 2 = iln(i? 2 /-Ri)- This case corresponds to the vortex chaine in the ring. 

If we have N pairs “vortex-antivortex” located uniformly on the circumference, one can interprete 
this as two vortex chaines. The complex-conjugated velocity has the form 


w'M = t~ 


c 



c 



2 iri r 

-m — 

U)\ R 2 


(6.5) 


where z\ (respectively zq) is the coordinate of a vortex (respectively antivortex) in the chaines 
and ni = 1, 77-2 = —1- For the sake of simlicity we can take z\ = ro and Z 2 = roexp(ni/N). 

To obtain the velocity of the jth vortex, we should substract the term corresponding to this 
vortex and then put z = Zj. We have 


W'(zj) = 


rik 


( I i In — 

Zk 


— ( ( i In 


ZjZ k 

Rl 


— ( ( 2 i In — 

S 1 R 2 


2ir> Tj 

-In — 

R,2 


( 6 . 6 ) 


where if j = 1, then k = 2 and if j = 2, then k = 1. 

Notice that the real semi-period of the (-functions is equal to uq = n/N while the imaginary 
semi-period does not change, i.e. LU 2 = i\n(R 2 /Ri). As in the part A we obtain that 


W’(zj) = 


77 1 + 77 2 


Z. 




]UJ2_ 

UJl 


77 1 + 77 2 


m 

2 c 


0, ?7 1 = -?7 2 , 

lN 

-, 77 1 = 77 2 = 1 


(6.7) 


for 7~i = r 2 = \fll\R 2 , Oj — 6 k = 7 t/N. This means that the stationary points for two chains 
of alternating-sign particles are 


Zk = r 0 e l7Tk/N , r 0 = n/Ai R 2 - k = 1,2,... ,2N, 
n k = (-l) fe+1 . 


It occurs that there exist the stationary points in the case when all the vortices have the same 
signs 771 = 772 = ■ • • = 77at = 1 and are on the same circumference r = ro- To find these points 
we have to solve the equation W'(zj) = 0, where W'(zj) is the velocity of j vortex obtained 
from (6.4) by substruction of the term corresponding to the vortex, 


W'(zj) = 


77 


l 


, , 0 •, Zj\ 2iri r 0 

-^ 2,b ‘TJ + ^ lr ‘T 2 


(6.9) 


We have 

/ \ °o 2 1 / \ 

coth ( N In — ) = 4 —- — ttt sinh ( 2tN In — ) , 

V RiJ h l ~ q V Rz) 

f 7TI Cl)*? I \ 

where q = exp-— is the denominator of the (-functions. 

V wi J 

This equation has a root for any N and any R\, R 2 . For example, if we take R± = 0.5, R 2 = 1.5 
and N = 2 (these parametres were used in the numerical simulations below) then ro = 0.651159. 
This value coinsides with the numerically calculated stationary point of the systems (^.4[), (|5.5|). 
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6.2 Motion ’’vortex-antivortex”: linear problems 


Let us investigate the motion of the pair “vortex-antivortex” in the ring, i.e. we have n\ = 
— n 2 = 1 and the stationary points are defined by (6.3). 

If the initial data are close to the stationary points one can describe the motion of two vortices 
in an analitic way. Taking into account that the stationary points correspond to the solution of 
the equation W'(zj) = 0 we obtain 


W\ Zj ) 


1 


D L 


d Fj o d Fj 0 d Fj 0 dFj 0 

m T (xi ‘ Xl) + aS (l2 “ X2) + ‘ 


where Fj is defined in ( |6.1|) , (x^yj), (x!j, y® ) are the stationary points, Z k = R2 exp(i0fc), and 


JX, 


9Fj 

dXn RlR2 


n k (e 2 - ei) + 2 Uje 3 + 


2 rjnj 


IT 


dFj 

d Vj 


XjTlk 

R\ R-2 


( e 2 - ei); 


dFj _ x k n k 

dxk R\R -2 

dFj _ x k n k 
dy k Il\ R 2 


2 rj 

ei + e 2 + — ) , j 


7T 


(e 2 - ei). 


Here we have taken into account that z± = x® = ro = \flk\R 2 . z 1 = x® = — ro = — \/-Rii? 2 , 
C'(z) = —V(z) [Q and 

■P In -jjjr ~ ~ 1 + w 2 ) = e 2 ; 


V (i In ^ - (0j - = V{uji) = ep 


P ( 2iln-^- ) = V(u ; 2 ) = e 3 


at the stationary points (x^yj), (x^y!^). Here ei, e 2 , e 3 are the roots of the characteristic 
equation (p]) and e\ + e 2 + e 3 = 0. Denoting 


a 1 = 


R1R2 

2 


a 2 = 


R 1 R 2 


(ei - e 2 ) 


2n" 

-e 3 + — ) , 

7T 


03 = 


R 1 R 2 


2 r] 


e -2 — ei - 2e 3 -, 

7T 


we can represent in the form 


W^(zi) = s [-01 (yi - y?) - ai(y 2 - y 2 ) - a 3 (xi - x°) + a 2 (x 2 - x°)] ; 


W'(z 2 ) = - [ai(yi - y°) + oi(y 2 - y°) - a 2 (xi - x?) + a 3 (x 2 - x°)] ; 


Denote x 3 — x® = x\, x 2 — x[] = x 2 ,... the deviations from the stationary points. Then we obtain 
the following systems 
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a) for the heat-flow equation 


' dx i 

—— = 03X1 — 02X2, 
at 

dx 2 
dt 

< 

dyi 
dt 

dy 2 _ . 

- 7 - = aiyi +aiy 2 , 
dt 

b) for the Schrodinger equation 


= —02X1 + CI3X2, 


— “iyi + 012/2, 


( 6 . 10 ) 


' dx 1 
dt 


— —a Wi — ai2/2, 


dx 2 

dt 


aiyi + aiy 2 , 


dyi 

dt 


— a 3 X 1 — 02 X 2 , 


#2 

—— = 02X1 — 03X2. 
dt 


( 6 . 11 ) 


Notice that e 3 < e 2 < 0 < e\ (f|), hence, 01 > 0, 02 > 0 and the sign of 03 depends on the 
parameters of the ring. We obtain the systems of ODEs with the constant coefficients. 

Let the initial positions of the particles be at the points ( 5i,£i ) (i = 1,2), where 5i,£i are small 
parameters. Then we can obtain an analitic form for the low of small motions near the stationary 
points. 


6.3 Heat-flow equation 


The system (|6.10|) splits in two independent systems for xi,X 2 and yi,y 2 and we have 


with 


Since 


xi(t) = C 1 e^ a ’ 3 ~ a2 ' )t + C 2 e^ +a ^\ 
x 2 (t) = - C 2 e^ 3+a ^ t 


C\ = 
C 2 = 


<5i + 82 


<5i — 82 
2 


03 + «2 


4 

R 1 B/2 


(e 2 - e 3 ) > 0 , 
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2 ( 4n\ 

03 ~ “ 2 = Rjh \ 2 ~ t ) < °’ 

we have e(“ 3+ “ 2 ^ —* oo and e^“ 3_a2 )* —> 0 as t —> oo. If the initial positions of the particles are 
symmetric with respect to the stationary points, <5i = —e >2 = <5, we obtain x\ (t) = —X 2 (t) = 
£ e (a 3 +a 2 )t p e _ particles move away from the stationary points. For antisymmetric initial 

positions 8 \ = 82 = 8 , and x\(t) = X 2 (t) = 8 e^ a3 ~ a2 ^ —> 0. This means that the particles 
described by the linearized equation are attracted to the stationary points. 

From the system for %){ we have: 

fyi(f) = £ 2 e 2ait , 

{V 2 (t) = £ 2 e 2ait - (£2 ~ £ 1 ). 

Both the solutions are divergent. If at the initial moment the particles are at the x-axis (ei = 
£2 = 0) they move along this axis. 

More detail numerical investigation of the motion of the pairs “vortex-antivortex” (not only for 
N = 2 but for any even N) shows that all the particles move from the stationary points to the 
nearest boundary or to the nearest neighbours and annihilate. For example, initial positions are 
far enough (the vortices are distributed uniformly on the same circumference), then they “don’t 
fill” each other, move to the boundary and annihilate with their reflections (see Fig.6a). 

If they are close enough, they attract and annihilate with the nearest neighbour (see Fig.6b). 
Such a behaviour is similar to the motion of the particles in the circle and corresponds to the 
usual physical picture: the charges with the different signs attract. 


Alternating signs, N=6 


Alternating signs, non-symmetric initial data 



Fig.6a. Even number of the alternative-sign Fig.6b. Particles are close enough, 

vortices 

If we have all the particles with the equal signs they always move to the outer boundary and 
annihilate. The pictire of motion is the same as for circle (see Fig. 5b). 
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6.4 Schrodinger equation 


In the case of the Schrodinger equation the system does not split and we have more complicated 
motion. We see that x 2 = x\ + C and the second part of the system takes a form 


dy 1 

at 


= a 3 xi - a 2 (C - xi), 


dm 

dt 


= a 2 x 1 - 03 ( 6 * - X\). 


Eliminating y\, y 2 from the equation for x\ we obtain 


d 2 x 1 


dt 2 


- 2 ai(a 2 + a 3 )xi - Cai(a 3 + a 2 ). 


In the standard way we find the solution of the system satisfying the initial data 5q(0) = <5i, 
x 2 ( 0 ) = S 2 , ?/i( 0 ) = £ 1 , and y 2 ( 0 ) = £ 2 : 


x\{t) 
x 2 (t) 


^ cos kt — d 2 sin kt + ^ , 

— ——— cos kt — d 2 sin kt + ^ . 


Here k equals k = y^2ai(a 2 + a 3 ) and d 2 is a constant that will be determined below. 
The system for y t with initialy conditions 2/1 (0) = E\ and 2 / 2 ( 0 ) = e 2 yields 

()= (*i-&)(tt2 + «3) in k £i±£l k (<?i + <? 2 ) (a 3 - «2> 
y w 2k 2 2 2 

»(t) = + - Mhi + °3) si „ a + £i±£i cos « + (fr + frK°2-°3) f _ £i_£i. 

We obtain the constant d 2 from the equation for x 2 : x 2 (0) = ai(ei +e 2 ). Namely, 

ai(ei + e 2 ) 


d 2 = — 


( 6 . 12 ) 


k 


We have periodic solutions for all the variables if d'i = — S 2 = 5 and £\ 
positions of the vortices are symmetric with respect to y-axis. Thus 

c 1 . ~ H a 2 + a 3 ) . 

x\ = 0 cos kt, 2/1 =---sm 

, r , ~ <^(02+03) . 

x 2 = —0 cos kt, 2/2 =- 7 -sm 


= e 2 = 0 , i.e. if the initial 

kt, 

kt. 


The vortices describe elliptic trajectories. 

In general any number of the pairs “vortex-antivortex” with symmetric initial positions moves 
along the closed trajectories (see Fig.7a for N = 2 and Fig.7b for N = 6 ). If the initial positions 
are close to the stationary points, the trajectories are ’’narrower” and ’’shorter”. 
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rl =(0.8,0.8), r2=(-0.8,0.8), n1=-n2=1 


Alternating signs, N=6 



Fig.7a. 


Fig.7b. 


Non-symmetric initial positions lead to a very complicated motion due to the terms which are 
linear in t (see (|6.12| )). The velocities of the vortices are different and the picture of the motion 
for the linear problem strongly depends on the initial data. 

Thus, if we have the particles with the different signs with arbitrary initial positions we obtain 
chaotic non-periodic motion (see Fig.8a for N = 5). Even the initial positions are symmetric 
but the number of particles is odd, the motion is chaotic (Fig.8b). 

Motion of 3 vortices: T=1.0 Motion of 5 particles: symmetric initial data 




Fig.8a. 


Fig.8b. 


Motion ” vortex-vortex” . The stationary points in this case is not so convenient to calculate 
the coefficients of the linear system. But it is obvious that if the initial positions are symmetric 
the vortices move by the same way. It means that the particles move independently in the same 
direction on the same circumference as it was for the circle (see. Fig.2b). If the initial positions 
are not symmetric, the vortices move globally in the same directions but along the comlicated 
trajectories with the loops (see, for example, Fig. 9 for N = 3). 
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Motion of 3 vortices: T=1.0 



Fig.9a. 


Motion of 3 vortices: T=2.5 



So, for Schrodinger equation we have very different pictures of motion depend on the initial 
positions of the vortices. 
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